Nonlinear coupling between excited atoms and surface polaritons 
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Casimir-Polder interactions between an atom and a macroscopic body are typically regarded as 
due to the exchange of virtual photons. This is strictly true only at zero temperature. At finite 
temperature, real-photon exchange can provide a significant contribution to the overall dispersion 
interaction. Here we describe a new two-photon resonant process between an atom and a planar 
interface. In order to understand recent experimental results with hot atomic vapours in glass 
microcells, we derive a quantum mechanical description in order to understand in detail the possible 
coupling between atoms and surface polaritons. We derive a nonlinear effective Hamiltonian to 
explain how atoms can couple resonantly to the surface polariton modes of the dielectric medium 
which leads to second-order energy exchanges. 

PACS numbers: 12.20.-m, 42.50.Nn, 71.36.-l-c 



I. INTRODUCTION 

Fluctuation-induced forces such as Casimir-Polder 
(CP) forces between atoms or molecules and macroscopic 
bodies are manifestations of the zero-point energy of the 
electromagnetic vacuum plj . They occur even if the atom 
and the macroscopic body are in their respective (unpo- 
larized) ground states [2] and can be understood — at 
least in the nonretarded Hmit — as interactions between 
a spontaneously generated atomic dipole and its (instan- 
taneous) image inside the macroscopic body. As soon 
as it became possible to achieve atom-surface distances 
below 100/im, experiments revealed that the couphng be- 
tween the atom and the surface at these short distances 
would produce significant effects [5- 

Following the advances in laser coohng and trapping 
techniques in the 1980s, a new area of research has 
emerged. Modern laser-based techniques have allowed 
an unprecedented amount of control, with this control 
came the ability to study very large atom-based systems. 
As a result of these advances, trapping and manipulat- 
ing single atoms, driving atoms into highly excited Ry- 
dberg states or creating Bose-Einstein condensates have 
become possible. New complex microstructures like atom 
chips allow one to trap, cool and manipulate ensembles 
of ultracold atoms in the vicinity of a surface [4]. 

Atoms and surface polaritons are very distinct quan- 
tum objects with different characteristics which make 
them suitable to perform different tasks. Atoms are very 
good candidates for storing and manipulating quantum 
information. The extremely promising results in the field 
of Rydberg atoms, both in ultracold atoms or in thermal 
vapours, shown that they make very good candidates to 
built quantum gates [6]. The renewed experimental 
interest in Rydberg atoms is due to the unique opportu- 
nities afforded by their exaggerated properties [7] which 
make them extremely sensitive to small-scale perturba- 
tions of their environment and to dispersion forces. Pre- 
vious work |5j showed that these properties includes mas- 



sive level shifts that a Rydberg atom experiences in close 
proximity of another atom or in the vicinity of a macro- 
scopic body, with shifts on the order of several GHz ex- 
pected at micrometer distances. 

Surface polaritons appear at the interface of two media. 
They represent particular solutions of the Maxwell equa- 
tions which correspond to waves propagating in parallel 
to the interface and whose amplitude decreases exponen- 
tially when moving away from the surface. They are ca- 
pable of interacting and be moved around on a surface, 
making them very attractive means of transporting quan- 
tum information from one point to another [S]. Upon 
taking advantage of the individual properties of atoms 
and surface polaritons and their different properties, it is 
possible to propose sophisticated quantum circuits |10J . 

Atom-polariton couplings lead to the (nonresonant) 
CP interaction between an atom and a planar interface. 
In the nonretarded limit, this interaction scales with 1/z^ 
[z is the atom-surface distance) [2 . Moreover, it has al- 
ready been shown that it is possible to turn the (usually 
attractive) CP interaction into a repulsive force by a res- 
onant coupling between a virtual emission of an atom and 
a virtual excitation of a surface polariton |llj . Similarly, 
it has been shown that the atom-surface coupling can 
drastically modify atomic branching ratios, because of 
surface-induced enhancement of a resonant decay chan- 
nel [12]. 

In this article we analyze a new type of near-field ef- 
fect involving surface polaritons inspired by the experi- 
ment of Kiibler et al. [13] with hot Rb vapour in glass 
cells. Their experiment indicated that a description of 
the atom-surface interactions should also include a non- 
linear coupling between the atomic transitions and sur- 
face polaritons. Their experimental results indicated that 
it should be possible for an atom to be coupled resonantly 
to the surface polariton modes of the dielectric material 
which leads to second-order energy exchanges with the 
atomic transition energy matching the difference in po- 
lariton energies. 
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The article is organised as follows. After briefly re- 
viewing the formalism of macroscopic QED in Sec. 
we derive an effective nonlinear atom-polariton coupling 
Hamiltonian is Sec. |III| and give concluding remarks in 
SecHVl 



II. BASIC EQUATIONS 

In the electric dipole approximation, the Hamiltonian 
that governs the dynamics of the coupled atom-field sys- 
tem can be written as [2] 

oo 

H= du! (fir huj f^r, Lo) ■ f{r, Lo) + ''^^hLUn Ann 


E(rA) . (1) 

It is expressed in terms of a set of bosonic variables 
f^(r, uj) and f (r, uj) that have the interpretation as ampli- 
tude operators for the elementary excitations of the sys- 
tem composed of the electromagnetic field and absorbing 
medium. They obey the commutation rules 

[fk{r,io)jl,ir',L,')]=6kk'Siuj-Lu')d{v-r'). (2) 

The relevant part of the Hamiltonian for our study is the 
interaction Hamiltonian 



^int — ~ ^nmd„m ' E(ryt), 



(3) 



in which the Anm denote the transition operators be- 
tween two dipole-coupled internal atomic states |m) 
and \n) with dipole transition matrix elements dnm — 
(n| d |m). The frequency components of the electric field 
operator 



(4) 



E(rA) = / dujEivA, uj) + h.c. 



are constructed via a source-quantity representation from 
the dynamical variables F(r,w) and f(r,a;) as 

E(r^, uj)^ J d\ Ge{vA,v, uj) ■ f (r, uj) . (5) 

The tensor Ge(ryi,r, oj) is related to the classical Green 
tensor G(rA,r,a;) by 



Ge(r,r',w) = i- 



-Im£(r',a;)G(r,r',a;). (6) 



The Green tensor itself is a solution of the Helmholtz 
equation 



V X V X — T£(r, 1^) 



together with the boundary condition G(r, r', a;) — ^ for 
|r — r'l — oo. The Green tensor obeys the useful integral 
relation 

/ d^s — Ime(s,a;) G(r,s,a;) • G*(s,r',iu) 

= ImG(r,r',a;) , (8) 

which follows directly from the Helmholtz equation ([t]) 
and which reflects the linear fluctuation-dissipation the- 
orem. 



III. EFFECTIVE ATOM-POLARITON 
COUPLING 

In this section, we derive the quantum mechanical de- 
scription for an effective nonlinear atom-polariton inter- 
action. The situation we envisage is depicted in Fig. [T] 
in which an atomic transition couples resonantly to two 
surface polariton modes of the dielectric material. This 
corresponds to second order energy exchanges with the 
atomic transition energy matching the difference in po- 
lariton energies. To illustrate our basic idea, we consider 



dielectric 
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FIG. 1: Scheme of resonance between the atomic transition 
and one surface polariton which leads to the creation of a 
second polariton f^i ~ cjio -I- 0,2- 

the interaction of an atomic transition of frequency Wiq 
between two eigenstates |0) and |I) with two surface po- 
laritons with corresponding center frequencies Oi and 0,2 
(r^i ^ O2) for whom the condition ili w wio -I- O2 is sat- 
isfled. The polariton resonance frequencies f2i and 0,2 
are assumed to be far from any other atomic transition 
frequency w,„„. 

Heisenberg's equations of motion for the dynamical 
variables and the atomic transition operators follow from 
the Hamiltonian ([T]) as 



f(r,Ci;,i) ~ — iwf(r,a;) 
yi^ e(r, w) A/j/dfc; • G* (ry!i, r, w) , (9) 



k,l 



G{r,r',uj) = S{r-r') (7) 



^ [AknAkm - imfcd„fc ) • ^{y a) ■ (10) 

k 
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Formal integration of Eq. ( 10 ) yields 

A (f) — p*"™"* A 



X Akn{t')dk,n-A^k{t')dnk\-'E(rA,t'), (11) 



which, inserted back into the equation of motion for the 
dynamical variables of the medium-assisted field, Eq. (|9| , 
yields the first iteration of the equations of motion for the 
dynamical variables as 



i{r,uj,t) ~ — «ajf(r,a;) + 



/lme(r,a;) 

V hTTer, 



J2 e'^"^'^' Arandmn ' G* {v A, V, Uj) 



1 C C C 2/2 

+ ^ E \ ^^'y ^^'y rfV^^Vlni£(r,c^)Ime(r',w')e*"-^*~*'^dmn-G*(r^,r,c^) 

X Akn{t')dkm - i„fe(i')d„fc] • [G(r^,r',c^') • f(r',t^',t') - G* (r^, r', c^') • i\r\u',t') 

I 



Equation (12 1 is now a nonlinear operator equation that 



is capable of describing resonant processes involving two 
polaritons. This is despite the fact that the original 
Hamiltonian ([T]) is bilinear in all operators. The effective 
nonlinearity appears as a consequence of the iteration. 
In order to pick out the resonant interactions from the 
equation of motion, we introduce as usual slowly vary- 
ing amplitude operators as f(r,cj,i) — f(r,w,i)e*"* and 

Amn{t) = Amn{t)e~^'^"'"* and apply the Markov approx- 
imation. This involves taking the slowly varying ampli- 
tude operators out of the integral at the upper time t. 



|1>- 
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FIG. 2: The transition from |1) to |0) is mediated by the 
virtual transitions via the state |fe). 



For our investigation one needs to keep in mind that 
the polariton spectrum is not continuous but consists of a 
quasidiscrete set of lines of midfrequencies Qk and widths 
7fe, where the linewidths are typically very much smaller 
than the line center separations jk ^ {^k+i — ^k-i)/2. 
We then divide the w axis into intervals ~ [{flk-i + 
ilk) /2, {^lk + ^k+l)/2]■ Recalling the resonance condition 
r2i ~ Wio + we apply the rotating- wave approxima- 
tion and finally arrive at the effective equation of motion 
describing the dynamics of the resonant atom-polariton 
coupling where now the frequency integrals have to be 
taken over the linewidth of the surface polaritons. 



Here we have introduced the abbreviation 



A 



10 \ ^ 

k 



g(r,r',w,a;') = 
dfeo ■ Gl{rA,r,uj)(g)dik ■ Ge(rA,r',a;') 



Ukl + uj' 

difc • G*(r^,r,w) (g)dfeo • Ge(rA, r', w') 



wofc + w' 



(14) 



for the operator-valued coupling tensor. As one can see 
from the structure of g(r, r', w, w'), the atom-polariton 
coupling is mediated by a virtual atomic transition from 
|1) — >■ |0) via an intermediate state \k) (see Fig. [2| with 
dipole moments dn: and dofc. The equation of motion 
( 12 1 can be thought of as being generated by the effective 



nonlinear interaction Hamiltonian 



f (r, bj) — — iwf (r, w) 



-I J (Fr' J dJ 

A,,/ 



(13) 



H^ff ^hj d^r' J d^r j duo j duj' 

A^, 

f^(r,a;) • g(r, r', w, w') •f(r',a;') - (lu ^ cu' , r ^ r')* 



(15) 
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where the second term describes the process that is ob- 
tained by interchanging uj with w', and r with r'. 



A. Coupling for singly excited polaritons 



The interaction energy (15 1 depends, by virtue of its 



dependence on the Green function, on the atom-surface 
distance z and thus contributes to the total Casimir- 
Polder potential. In the following, we focus on the ef- 
fective potential that results resulting from this atom- 
surface coupling. For this purpose, it is useful to define 
single-polariton excitations (similar to the construction 
in Ref. [14]) as 



|l(rA,f^.)) = 



xG:(r^,s,c.).ft(s,c.)|{0}) 
with the normalization factor 



(16) 



the form 



C/cff 



71 72 



Tr ImG(r^ , r ^ , i ) Tr ImG(r ^ , r a , 2 ) 
^ |Tr[ImG(rA,rA,f^i) • dofe (g)dfci • ImG(rA, r^, ^^2)] 



(f2i+c^ofe)2+7iV4 
-Tr [ImG(rA, r^, ^1) ■ d^i «) dofc • lmG{rA, r^, ^2 



(20) 



Comparing Eq. ( 20 ) with the Casimir-Polder potential 



at finite temperature [H |H] , 



5(i"A,f^i) 



Mo 



r!|TrImG(rA,rA,r!.). 



(17) 



Using the integral relation ( |8| f or the Green tensor, ine 
easily checks that the states ( |16| ) are indeed properly nor- 
malized. 

Note that the states |l(rA,r2i)) carry a vector index 
as well as the continuous space and frequency labels. 
In this approach we assume that the local density of 
states w^Im G(r, r, w) can be approximated by a single 
Lorentzian peak of mid-frequency Q and width 7, 



cj^Im G(r, r, lu) ~ f^^Im G(r, r, n) 



7V4 



(w- 17)2 -f 72/4 



(18) 

In our envisaged situation of a resonant coupling between 
a single atomic transition and the difference between two 
polariton resonances, the energies of the initial and final 
states are identical. Degenerate first-order perturbation 
theory asserts that the interacting potential is [15] 



C/cff 



{K\ \I) 



(19) 



Here |/) — \Ia) \0i) II2) stands for the tensor product of 
the initial excited atomic state \1a) and a singly excited 
polariton with frequency V,2, and \K) = \0a} |1i) IO2) de- 
notes the tensor product of the final atomic state jO^) 
and a single excitation in the polariton with frequency 
ill. The single-polariton states |lj,) = |l(rA, fi^)) are de- 
fined according to Eq. ( [l6| and |0,y) denotes the polariton 
ground state \0^) = \{0Y), yuj e [fl^ - 6uj/2, + Su}/2]. 

Using the commutation rules of the operators as well as 
the properties of the Green functions, together with the 
definition of the Lorentzian lineshape, Eq. (18), we find 



that the effective interaction potential can be written in 



AijNR = ^ofcsT^' ejTr [a(*0) • G(rA, r^, i^,)] 
3=0 

+/io X! w^„nth(wfc„)d„fc • Re G{rA, r^, Wfcn) • dfc„. 



(21) 



where are the Matsubara frequencies, one observes 
that the effective potential dependence on the atom- 
surface distance z is given by the same (inverse cubic) 
power-law as the nonresonant Casimir-Polder potential. 
This is because the second-order nonlinear Hamiltonian 



( 15 1 enters only a (degenerate) first-order perturbation 
theory (19), whereas the nonresonant potential arises 



from a bilinear Hamiltonian in second-order perturbation 
theory [2^. 

The total potential experienced by the atom is the sum 
of the nonresonant (attractive) CP potential and the res- 
onant coupling between the atoms and the surface po- 
laritons, 



AE 



iTotal 



AE 



NR 



AE' 



(22) 



The respective energy shifts for the nonresonant and 
(second-order) resonant interactions in the nonretarded 
limit are 



UttHz^ 



-^^^nth(t^te)|d„fc|'Re 



e('^fcn 



.^(W/c„) 



(23) 



5 



and 



71 72 



2z3 Y TrImG'(f7i)TrImG'(r22) 
^ jxr [ImG'(r!i) • dofc ® dki ■ ImG'(r!2)] 



ill + t^Ofc 



(fli +cjofc)^+7iV4 
-Tr [ImG'irii) ■ d^i (g) do^ • ImG'(f72) 

''(lli+c.fei)'+7iV4j' 
where we have defined 

1 

,0 2; 



327rl72 



£(^]) - 1 

e(fl) + 1 



(24) 



(25) 



as the nonretarded scattering Green tensor of a planar 
halfspace in the coincidence hmit r ^ r', G'(fi) = 
z^G{r,r,n). 



B. Thermal States 

As we are deahng with thermaUy excited surface po- 
laritons the concept of perturbation theory has to be ex- 
panded from pure states described by a single state vector 
to a statistical mixture or ensemble of states [T5]. The 
density matrix for a thermal state with temperature T 
can be written in the Fock basis In) as 



Pth =^Pn 



where Z{T) = 



E 



-nhn„/kBT 

Z{T) 



(26) 



-mtinr^/kBT (jgnotes the partition 



function. 

So far we have computed the interaction energy for 
the situation in which there is initially only one excited 
polariton with frequency 0.2 and in the final state only 
one polariton with frequency ili [see Eq. (19)]. This has 
to be generalized to thermal states in which there can be 
initially m polaritons with frequency and n polaritons 
with frequency $72- In this case, we rewrite the result of 
the perturbation theory as 



{K\H,s\I) =Tr H,sknH.s\K) {K\ 



(27) 



where 



Pin = PthiP-l) ® Pth{^2) ® , 

= E?'™ ' \mun2) (mi,n2| ® \Ia) {Ia\ ■ (28) 

Due to the form of the effective interaction Hamiltonian 
iJcff the only final state \K) that provides a non- vanishing 



matrix element will be \K) = |(m+ l)i. (n — 1)2) |Oyi)- 
Hence, 



{K\H,s\I) = ^pWp(fn"^ + l)(i)(")(2)C/e'ff, 



= [nth(f^l) + l]nth(f^2)C/cff 



(29) 



Finally, the resonant energy shift AE^ for thermal states 
will be given as 



AE'- 



Ucs\/[nthifli) + 1] nth(fl2 



(30) 



This result is intuitively clear, as the initial polariton 
with frequency ^2 has to be thermally populated before 
the resonant interaction can take place. 



IV. FINAL REMARKS 

The key idea in current experiments is to confine ther- 
mal vapours of — typically alkali — atoms in spec- 
troscopy microcells whose thickness lies in the 20 — 
1000 nm range. These vapour cells are typically made 
of (transparent) dielectric materials such as sapphire, 
quartz or fluoride crystals. The nonretarded nonlin- 
ear interaction potential between an atom and a dielec- 
tric surface is governed by the (nonnegative) dissipative 
part Imrp(aj) of the surface response function, where 
rp{u!) = (£(aj) — l)/(e(a;)-|-l) is the Fresnel reflection coef- 
ficient for p-polarized waves at a planar interface between 
a dielectric with permittivity £(cj) and free space [see 
Eq. (25)]. For each dielectric material there is a unique 



surface polariton spectrum that depends sensitively on 
the concentration and distribution of the impurities and 
surface quality of the samples (surface roughness). As 
each sample is unique, the polariton spectrum should be 
found experimentally. The Fresnel coefficient rp{uj) ex- 
hibits resonances (strictly speaking, there are poles in the 
complex frequency plane where e{u!) = —1), the associ- 
ated modes are the surface polaritons that are linked to 
the bulk resonances of the dielectric permittivity e(w). 
We use a multi-resonance Drude-Lorentz model 



e{uj) = 1 + 



''Pi 



(31) 



with plasma frequencies ujpj and transverse resonance 
frequencies ujxj as an approximate description of the di- 
electric permittivity. Selecting two well separated reso- 
nances, the effective potential becomes 



AE 



R 



flQC UJpiUJp2 



1287rz3 ^/njh 

r— — , . , _ — 77;r^v^ 5 [dofcl |dfci| 
X V[?^th(^ll) + 1] fith(fl2) 2^ 

k 



12 



r^i + wofc 



(J7i+LJofc)2 + 7iV4 (fli +Wfei)2 +72/4 



(32) 
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Closer inspection reveals that the interaction is typi- 
cally dominated by the nonresonant CP potential, whose 
maximum value for a one-polariton model is 



AE- 



■NR 



Hoc^ ksT ujp \dii, 



487rz3 h 



E 



'247rz3 



^nth{uj^i) \di^\ 



xRe 



(33) 



Note that we have considered an isotropic atom in this 
approach. 

Finally, let us investigate which intermediate atomic 
transitions might provide the largest nonlinear effect. In 
order to have an effect that is comparable to then non- 
retarded Casimir-Polder interaction, there has to be a 
matching atomic transition between energetically close 
states — note that the transition |I) — >■ |0) does not 
need to be allowed by the selection rules — i.e. the in- 
termediate state \k) has to be close to the initial and 
final states |1) and |0). The reason for this constraint is 
the rapidly decreasing magnitude of the dipole transition 
matrix element between states with increasing energy dif- 
ference. Therefore, the dominant contribution will come 
from an intermediate state | k) approximately halfway be- 
tween initial and final states. 

In this case the difference between the resonant and 



nonresonant terms will come from the last line in 
Eq. (32 1. Its maximum value is obtained whenever wofc 



or Wife is one of Oi ± 7i/2; away from these points the 
numerical value of this term decreases. As we have as- 
sumed throughout our calculations that all atomic tran- 
sitions are far from any single-polariton resonance, the 
Lorentzian peaks have to be broad, i.e. 71 has to be 
large. This in turn means that, in order for this non- 
linear effect to be comparable to the nonresonant CP 
potential, a strongly dissipative material is needed. For 
typical cell materials such as sapphire jl7| and quartz |18| 
the resonant nonlinear shift is too small to be relevant for 
current experimental investigations (only several kHz for 
Rydberg atoms). In order to access this nonlinear effect, 
future work will have to focus on investigating suitable 
cell materials for specific atoms. 

Although cell materials that are used in current exper- 
iments have usually been studied in depth, experimental 
findings present large discrepancies, due to variations in 
the quality of the sample, the degree of its impurities 
and the oritentation of the crystal axes [19] . Proper (ex- 
perimental) characterization of the cell material and its 
resonances is a crucial step in modeling this process. 
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